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Abstract - This paper presents an approach for reducing 
speckle in ultrasound images using Coupled Partial 
Differential Equation (CPDE) which has been obtained by 
uniting second-order and the fourth-order partial 
differential equations. Using PDE to reduce the speckle is 
the noise-smoothing methods which is getting attention 
widely, because PDE can keep the edge well when it 
reduces the noise. We also introduced a median regulator 
to guide energy source to boost the features in the image 
and regularize the diffusion. The proposed method is 
tested in both simulated and real medical ultrasound 
images. The proposed method is compared with SRAD, 
Perona Malik diffusion and Non linear coherent diffusion 
methods, our method gives better result in terms of CNR, 
SSIM and FOM. 

Index Terms - Ultrasound, Speckle, PDE, median filter, 

I. Introduction 

ULTRASOUND speckle is the result of the diffuse 
scattering, which occurs when an ultrasound pulse 
randomly interferes with the small particles or objects 
on a scale comparable to the sound wavelength. 
Speckle is an inherent property of an ultrasound image, 
and is modeled as spatial correlated multiplicative 
noise. In most cases, it is considered a contaminating 
factor that severely degrades image quality. To improve 
clinical diagnosis, speckle reduction is generally used 
for two applications: visualization enhancement and 
auto-segmentation improvement. Most speckle filters 
are developed for enhancing visualization of speckle 
images [l]-[3] including temporal averaging, median 
filtering, adaptive speckle reduction (ASR) and 
wavelet based methods [4] -[5], frequency 
compounding methods [6]- [7]. 

However, above methods could not succeed to 
balance between speckle suppression and feature 
preservation due to the complexity of speckle statistics. 
The comparative analysis [8] and review methods [9] 
also says that a technique that relies on a more accurate 
model for removing speckle noise while preserving 
image feature well would be rather valuable for 
practical use improve image simplification, which is in 
turn very beneficial in automated object detection 

II. BACKGROUND OF PDE 

The use of partial differential equations (PDE) in 
image processing has grown significantly over the past 
years. Its basic idea is to deform an image, a curve or a 



surface in a partial differential equation framework, and 
to approach the expected result as a solution to this 
equation. 

P. Perona and J. Malik put forward an anisotropic 
diffusion (AD) equation to smooth a noisy image: 



du(x ,y , t) 
It = 



div [g\Vu (x,y,t)\)Vu(x,y,t) 



(i) 



where du(x,y,t) :£2,X— >(0,+oo)— ►.R is a scale image, g( | 
Fu | ) is a decreasing function of the gradient. Their 
work made an important influence on this field. 

2D SRAD takes the format of the PDE of 
conventional anisotropic diffusion. Given an intensity 
image Io(x, y) having none zero-valued intensities over 
the image domain U, the continuous form of SRAD is 
expressed as follows [10] 



h 



div[c{xj,:t)U(x,y,t)] 



7(xj:0)=/ (xj),(J/(xj;0/Jfl = 



(2) 
where d£l denotes the border of Q, n^ is the outer 
normal to the dQ, c(x) is the diffusion coefficient, and q 
is the ICOV. 

Both P-M diffusion and SRAD introduce blocky 
effects in image. So Yu-Li and M.Kaveh [11] proposed 
a fourth-order partial differential equation. This 
equation can reduce the noise, at the same time it can 
keep the image edge better. The PDE is: 



<LL. 

dt~' 



-V 2 [s(|V 2 /| 



Where g\ 



1 



-V[g(*)] 

f~ , k is constant. 



(3) 



[l + (slk) _ 

The advantage of using fourth-order PDE in denoising 
image is it removes the blocky effects that made by 
second-order nonlinearity diffusion equation. 

III. PROPOSED METHOD 

The coupled PDE based on second- order and fourth- 
order PDE is : 

d -±=-V[g x {AI)V I + a.g 2 {AI)AV I] (4) 
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Where 'a' is the coupled coefficient and the functions 
gi, g2 are plus and non-increasing function 

The median filter is a well-known "edge preserving" 
nonlinear filter. It removes the extreme data while 
producing a smoothed output. We propose to use a 
median regulator in coupled PDE, which combines 
both medina filtering and strength of unity of second 
order and fourth order Partial Differential equation. 

The proposed technique is defined by the following 
relations: 

dl(x,y,t) 



dt 



■ = -V[ gl (AI)VI + 

a.g 2 (AI)AVI] + af 



(5) 



l{x,y,t) t=0 = I 
d n l—0,f—median(l). 

Speckle noise is signal-dependent noise. Typically, 
the bright regions have stronger noise than the dark 
regions. The source term /provides a boosting force to 
guide (or normalize) the diffusion evolution and also 
accelerate the convergence rate compared to natural 
diffusion, a is control coefficient whose value lies between 
Oand 1 

The function gi, g2are defined as 



g l {x) = exp{-[x/k l ] 2 ), 
1 



g 2 ( x ) = 



(6) 



(l + (x/£ 2 ) 2 ) 



Where ki, k 2 are constants. The coupled PDE removes 
the noise and keeps the edges, regulator act as an 
accelerating force to produce an image with less noise 
and enhanced structure 

IV. IMPLEMENTATION AND RESULT 

The proposed method is tested on a simulated fetus 
image using Field-II ultrasound simulation software 
and a real ultrasound liver image. 

Throughout this paper, we take multiplicative noise 
as the model for speckle. The degraded intensity for 
each pixel will be 

/. .— f. u. . (7) 

i, j J i,] i,j \') 

Where f t - is the intensity at pixel (i,j )in the 
original image, U { . the multiplicative noise, and 
1 1 • is the intensity at pixel (i , j ) in the observed 
image 

The structure similarity (SSIM), Contrast to noise 
ratio (CNR), Pratt's Figure of Merit (FOM) are 
calculated to compare the performance of proposed 
method with SRAD, Perona Malik diffusion and Non 
linear coherent diffusion methods. 

The structure similarity (SSIM) [12] was used to 
evaluate the overall processing quality 



SSIM = - 



(2^^+C 1 )(2a xj; + C 2 



-// 2 +C 1 



)(* 



2 + <7 2 +CJ 

y 2 ; 



(8) 



[My. ' ' ^ l 

Where fi t and o t are the mean and standard 
deviation of intensities of pixels in a local window, and 
is the constant to avoid instability. All involved 
parameters were set as suggested in [13]. The value of 
SSIM varies between andl, with unity for ideal 
processing quality, i.e., when the filtered image is equal 
to the reference image. 

The contrast-to-noise-ratio (CNR) or lesion signal- 
to-noise ratio [13], [14], was computed by 
\ju,-ju 2 \ 
CNR =7= (9) 

2 

Where ju l and a l are the mean and variance of 

intensities of pixels in a region of interest (ROI), and 

2 
1^2 and a 2 are the mean and variance of intensities 

of pixels in a background region that has the same size 
as the ROI to be compared with. 

Edge preservation or edge enhancement can be 
quantitatively evaluated by Pratt's figure of merit 
(FOM). 

N 

I 



FOM=- 



1 



1 



max \N ,N 



ideal 



i=l \ + d.k 



(10) 



Where TV, N ideal are the numbers of detected and 

ideal edge pixels, respectively, di is the Euclidean 
distance between the /th detected edge pixel and the 
nearest ideal edge pixel. A is a constant typically set to 
1/9. The dynamic range of FOM is between [0, 1]. 
Higher value indicates better edge matching between 
processed image and the ideal image. 

In the implementation of median regularized coupled 
PDE method the following values are selected for the 
parameters: ki=0.03, k 2 =0.04 and a=l.The mask size of 
3X3 is chosen to calculate median regulator. 

The result of proposed method, NCD, SRAD and 
Perona Malik diffusion on Field II ultrasound 
simulation software generated fetus image are shown in 
fig. 1. The same methods are tested on and real 
ultrasound liver images and the results are shown in 
fig.2.Comparision of performance measure are shown 
in Table. 1. Image profile along 154 th row of various 
filtered images are shown in the fig. 3 



V. 



DISCUSSION AND CONCLUSION 



The result processed by the new method is much 
sharper in terms of edge preservation and smoother in 
terms of speckle noise reduction than the other three 
filtered results. The execution time is also much shorter 
than the other methods. For quantitative quality 
evaluation, we provide three metrics SSIM, CNR and 
FOM. Table. 1. shows the evaluation results for the 
processed image. 
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The FOM value indicates that the new method is 
better than other methods in terms of edge preserving 
ability. CNR and SSIM values indicate that the new 
method gives a better processing result in terms of 
structure preservation and contrast enhancement. 

In order to keep the image feature better we adopted 
coupled PDE by uniting second order and fourth order 



PDE, Further more we added a median regulator to 
accelerate the convergence rate. The proposed Median 
regularized coupled PDE provides improved result over 
other methods in terms many parameters. Therefore it 
has a wide applicability in practice. 




Fig.l Result on simulated image (a) Original Noisy Image , Image filtered by (b) Perona Malik Diffusion 
(c) Non Linear Coherent Diffusion (NCD) (d) SRAD (e) Medial Regularized Coupled PDE 

Table 1 Comparison of Performance Measure 





NOISY IMAGE 


Perona Malik 


NCD 


SRAD 


Proposed 


CNR 


0.52±0.07 


0.88±0.47 


1.52±0.007 


1.77±0.37 


2.56±0.46 


SSIM 


0.4322±0.0028 


0.7279±0.0031 


0.7422±0.0035 


0.7531±0.0039 


0.8122±0.0038 


FOM 


- 


0.6543 


0.5873 


0.6546 


0.9067 




Fig.2. Result on real ultrasound image (a) Original Noisy Image , Image filtered by (b) Perona Malik Diffusion 
(c) Non Linear Coherent Diffusion (NCD) (d) SRAD (e) Medial Regularized Coupled PDE 
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Fig.3 Profile of the filtered image along 154 th row (a) Perona Malik Diffusion (b) Non Linear Coherent Diffusion (NCD) (c) SRAD 

(d) Proposed Medial Regularized Coupled PDE 
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